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ABSTRACT 

We study a large sample of narrow-line radio galaxies (NLRGs) with extended 
radio structures. Using 1.4 GHz radio luminosities, L1.4, narrow optical emission 
line luminosities, L[oni] and Ln a , as well as black hole masses Mbh derived 
from stellar velocity dispersions measured from the optical spectra obtained with 
the Sloan Digital Sky Survey, we find that: (i) NLRGs cover about 4 decades 
of the Eddington ratio, A = L bo i/L Edd oc L Une /M B H] (ii) L 1A /M BH strongly 
correlates with A; (iii) radio-loudness, 1Z = L 1A /Li ine , strongly anti-correlates 
with A. A very broad range of the Eddington ratio indicates that the parent 
population of NLRGs includes both radio-loud quasars (RLQs) and broad-line 
radio galaxies (BLRGs). The correlations they obey and their high jet production 
efficiencies favor a jet production model which involves the so-called 'magnetically 
choked' accretion scenario. In this model, production of the jet is dominated 
by the Blandford-Znajek mechanism, and the magnetic fields in the vicinity of 
the central black hole are confined by the ram pressure of the accretion flow. 
Since large net magnetic flux accumulated in central regions of the accretion 
flow required by the model can take place only via geometrically thick accretion, 
we speculate that the massive, 'cold' accretion events associated with luminous 
emission-line AGN can be accompanied by an efficient jet production only if 
preceded by a hot, very sub-Eddington accretion phase. 
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1. INTRODUCTION 

It became clear already shortly after the discovery of first radio galaxies that their 
strong radio emission is associated with a presence of luminous optical emission lines (Bade 
& Minkowski 1954; Osterbrock 1977; Grandi & Osterbrock 1978). Such an association was 
later confirmed by finding a correlation of radio luminosities with narrow line luminosities 
in Fanaroff-Riley type II (FRII, Fanaroff & Riley 1974) radio sources (Baum & Heckman 
1989; Saunders et al. 1989; Rawlings et al. 1989; Rawlings & Saunders 1991; Zirbel & Baum 
1995; Willott et al. 1999; Buttiglione et al. 2010; Koziel-Wierzbowska & Stasihska 2011: 
KS11). Using the narrow-line luminosity as a proxy for the cold accretion disk luminosity 
and the radio luminosity as a proxy of the jet power, Rawlings & Saunders (1991) found an 
approximate proportionality between these two quantities, with jet powers approaching and 
in some cases, even exceeding the bolometric luminosities of the accretion disks. 

Studies of the correlation of the jet powers with properties of central engines using radio 
and optical luminosities became more thorough and robust once methods of the black hole 
mass estimations have been developed (see, e.g., Woo & Urry 2002 and refs. therein). This 
allowed a determination of the properties of radio-loud AGN as a function of their Edding- 
ton ratio, A, defined as the ratio of the accretion bolometric luminosity to the Eddington 
luminosity. In particular, Koziel-Wierzbowska & Stasihska (2011) using the black hole mass 
estimates from the stellar velocity disperssion (a*) - BH mass (Mbh) relation (Tremaine et 
al. 2002) found that FRII radio galaxies span about 4 decades of the Eddington ratio and 
that, when scaled by their black hole masses, their radio luminosities correlate with Edding- 
ton ratio, in similarity to the correlation of absolute radio luminosities vs. absolute narrow 
emission line luminosities. One might consider this similarity as resulting from the fact that 
distribution of the black hole masses of FR II sources is rather narrow, with majority of them 
within the range 10 8 — 1O 9 M . However, since the physics of the jet production is linked 
to the Eddington-scaled accretion rate rather than to its absolute value and powers of jets 
produced with the same efficiency scale with the black hole mass, the 'primary correlation' 
to be considered should be the one between the Eddington-scaled luminosities. 

Such a correlation, however, contradicts the predictions of jet models which relate the 
strength of the central poloidal magnetic fields (those threading innermost portions of accre- 
tion disk and the black hole) with maximal pressure in the disk. According to the standard 
accretion disk theory, innermost portions of disks accreting at a rate corresponding to the 
Eddington ratio in the range 10~ 4 — 1 are radiation pressure dominated. Since this pressure 
does not depend on the accretion rate, the powers of jets - believed to be initially dominated 
by Poynting flux - are not expected to depend on the accretion rate either. In consequence, 
such models predict maximal jet powers ~ 100 times smaller than those observed in RLQs 
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with extended radio structures (see, e.g., Ghosh & Abramowicz 1997). 

The model which can account for the energetics of the most powerful jets and explain 
the observed radio - optical correlation is the one based on the so-called magnetically ar- 
rested/choked accretion flows (Narayan et al. 2003; Igumenshchev 2008; McKinney et al. 
2012). In such a model, the amount of the net magnetic flux amassed in the central region 
is so large that innermost portions of accretion disks are dynamically affected by central 
magnetic fields. In this regime, the accretion onto a black hole proceeds via interchange in- 
stabilities (= magnetic Ray leigh- Taylor instability: Stone & Gardiner 2007 and refs. therein) 
and magnetic flux threading the black hole is supported by the ram pressure of the accreting 
plasma. 

In order to directly confront this model with observations, we expand the sample of 
radio galaxies studied by KS11 by including other types of radio morphologies, investigate 
the correlation of radio- loudness 1Z (radio-to-optical luminosity ratio) vs. A, and investigate 
the dependence of source sizes on the Eddington ratio. The paper is organized as follows: 
in Section 2 we describe our sample selection and data reduction and analysis; in Section 
3 we present results of our analysis of optical and radio based correlations; in Section 4 
we investigate a consistency of these results with the jet production model which involves 
'magnetically-chocked' accretion scenario. Our main conclusions are listed in Section 5. 

Throughout the paper we assume a ACDM cosmology with H = 71 km s _1 Mpc" 1 , 
Q m = 0.27, and Q A = 0.73. 

2. THE DATA 
2.1. The sample 

To select our sample, we proceeded in a manner analogous to that described in KS11. 
Since we are primarily interested in radio galaxies with elongated structures, we started with 
Cambridge radio catalogs accessible by Vizier where such objects are well defined: 3C (Edge 
et al. 1959; Bennett 1962), 4C (Pilkington & Scott 1965; Gower et al. 1967), 5C (Pearson 
1975; Pearson & Kus 1978; Benn et al. 1982; Benn & Kenderdine 1991; Benn 1995), 6C 
(Baldwin et al. 1985; Hales et al. 1988, 1990, 1991, 1993ab), 7C (Hales et al 2007), 8C (Rees 
1990; Hales et al. 1995) and 9C (Waldram et al. 2003). We considered all the radio sources 
from these catalogs and cross-identified them in an automatic fashion with the sample of 
926246 galaxies with optical spectra from the Sloan Digital Sky Survey (SDSS) DR7 main 
galaxy sample (Abazajian et al. 2009). Taking into account the sometimes large positional 
uncertainties in the Cambridge radio catalogs, we adopted a maximum distance between the 
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position given for the radio source and that given for the optical galaxy varying from 0.2 
arcmin to 1 arcmin, depending on the catalog. We thus obtained a list of 2633 radio galaxies 
with available SDSS spectra. The final identification was done using the NVSS (Condon et al. 
1998) and FIRST (Becker et al. 1995) radio maps, which have much better spatial resolution 
than the Cambridge catalogs, and have been obtained at the same frequency of 1.4 GHz. 
For all the pre-selected objects we constructed jpg images superimposing NVSS or FIRST 
contours on SDSS images in the r band. The images are centered on the galaxy which is 
supposed to correspond to the radio source. By visual inspection of all these images, we found 
that 307 cases were actually misidentifications (among which many of them corresponding 
indeed to a galaxy but for which the SDSS spectrum was that of another, nearby galaxy), 
and 14 corresponded to spiral galaxies whose radio emission is produced over the entire 
disk, most likely by star formation. For the remaining 2042 radio sources we carried out a 
morphological classification by eye resulting in the following subclasses: 

• FR I type radio sources, considered as those where the maximum brightness of the 
lobes is closer to the center than to the extremity, 

• FR II type, where the maximum brightness of the lobes is closer to the extremity, 

• FR I/II type, where one lobe is of FR I type and the other of FR II type, 

• double-double radio sources where two pairs of coaxial lobes are detected, 

• X-shape radio sources with two pairs of lobes forming an X-shape structure, 

• one-sided radio sources showing only one lobe, 

• "elongated", i.e. radio sources that do not appear point-like, but whose angular size 
does not allow us to classify them more accurately, 

• radio sources, for which no morphological class could be assigned because of an atypical, 
irregular shape 

• compact, unresolved radio sources. 

This morphological classification, which is purely subjective, was carried out independently 
by two of the co-authors, DKW and GS. For about 15% of the objects we considered, our 
initial classifications diverged, although we agreed on our final classification. For further 
considerations in this paper, we restricted our sample to objects that clearly show the pres- 
ence of radio structures associated with radio lobes and/or jets, i.e. FRII, FRI, FRII/FRI, 
double-double, X-shape and one-sided. 
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We excluded the radio sources whose parent galaxies have a redshift larger than 0.4, 
in order to insure that the Ha emission line - which is crucial for our study - falls within 
the SDSS spectral range. We also excluded those radio sources where a broad component 
was clearly seen in the hydrogen emission lines. This allows us to make more accurate 
computations of narrow line luminosities as it avoids issues related to the decomposition of 
the narrow and the broad components. In addition, this limits our sample to objects which 
according to the Unified Scheme are observed at large inclination angles and by this step, we 
minimize viewing angle biases in our sample, in particular regarding source sizes. Finally, 
objects corresponding to galaxies where no emission lines were detected (after processing 
with the STARLIGHT code, Cid Fernandes et al. 2005, see below) were obviously removed 
from the sample, since all the considerations in this paper make use of emission line fluxes. 
After all these cuts, the entire sample on which this paper is based consists of 404 objects. 
They are listed in Table 1 (accessible only on-line), together with the properties that will be 
used. 

There are some differences between our starting sample and the catalog of elongated 
radio sources published by Lin et al. (2010). The latter was assembled via cross correlation 
of the SDSS DR6 with NVSS and FIRST. Those radio surveys have the advantage of being 
more homogeneous and deeper than the Cambridge catalogs. Indeed, the limiting radio flux 
density at 1.4 GHz in Lin's et al. sample is 3 mJy while that of the 3C catalog (when rescaled 
to 1.4 GHz) is 2 Jy, and those of the 4C-9C catalogs range between 400 mJy and 20 mJy 
(except for the 5C catalog which reaches 1.5 mJy but in very limited zones of the sky). One 
may wonder how many objects we are missing by using the Cambridge catalogs as a starting 
point rather than the NVSS one. It turns out that the number should not be very large, 
since our sample is limited to a redshift of 0.4, contains only objects with FRI, FRII and 
related morphologies which are the most luminous radio sources and, in addition, includes 
only objects with emission lines. All our objects have 1.4 GHz luminosities larger than 10 24 
W Hz -1 . In Lin's et al. catalog, there are only 120 objects with radio luminosities smaller 
than that out of a total of 1040 objects. This implies that, in our sample, we miss at the 
very most 12% of objects on the low luminosity side. On the other hand, for some reason, 
our starting sample contains 112 objects that fulfil all the selection criteria of Lin et al. but 
do not appear in their catalog (such as, e.g. the well-known objects 3C198, 4C+00.56 which 
are bright in radio and have large angular sizes). Another difference between our approach 
and the one of Lin et al. is that, while they attempted to define an objective way to trace 
the galaxy population smoothly from FRI sources to FRII, (as opposed to a sharp and 
perhaps arbitrary distinction between type I versus type II), we adhered to more a classical 
morphological classification paying attention to other types than just FRI and FRII. 
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2.2. Data processing 

The 1.4 GHz radio luminosities, L1.4, were obtained for each source from the sum of 
1.4 GHz fluxes of each of its components listed in the NVSS catalog, including the central, 
compact point source. The NVSS catalog was used here in order to avoid a flux loss from 
the extended and faint components missed in the FIRST catalog. 

The determination of the angular sizes of the radio sources depended on the morpholo- 
gies. For FRII radio sources, the angular sizes were defined as the distances between the 
hot spots or the most distant bright structures in opposite lobes. They were estimated from 
the FIRST maps, if available, or from the NVSS maps otherwise. The determination of the 
sizes of FRI radio sources is less straightforward, since lobes and plumes of the FRI radio 
galaxies fade away with the distance from the radio core. For each FRI source in our sam- 
ple, we obtained the 3 rms contour using the FIRST map and we took for the source size 
the largest extent of this contour, measured along a straight line crossing the radio core (if 
present). This procedure works well with straight sources. In the case of bent sources, the 
size determined in that manner is an underestimate. For one-sided sources, we adopted a 
procedure similar to the one for FRI sources to determine the angular size of the only lobe. 
The linear (projected) sizes were then determined from the observed angular sizes using the 
redshift of the corresponding galaxy as obtained from the SDSS. To meaningfully compare 
the sizes of all the sources, we divided by two the sizes of all the two-sided radio- sources: 
the result is called the characteristic lobe-size in the remaining of this work. 

The line luminosities of the galaxies associated with the radio sources are fundamental 
for our work. It is therefore important to determine them in the best possible way. For a 
removal of the stellar features from the observed optical spectrum, the best approach is to 
fit the observed continuum - excluding the spectral zones where emission lines are expected 
- with a composite stellar population obtained by spectral synthesis and subtract it from 
the entire observed spectrum. What remains is the pure emission line spectrum, whose 
intensities can then be measured. As in KS11, we have taken the [O m]5007, Ha and H/3 
line fluxes from the STARLIGHT databasd^, where they have been obtained precisely in this 
manner. 

As argued in KS11, we consider Ln a to be a much better measure of the bolometric 
luminosity of the AGN than the commonly used Inoni], because it does not depend on the 
ionization state. This has also been shown by Netzer (2009). The plots presented below are, 
however, presented in pairs, with one plot using L^ a and the other using L[oiii], for easy 
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comparison with works by other authors. 

In only about 20% of objects in our sample are both the Ha and H/3 line fluxes measured 
with sufficient accuracy to allow a meaningful estimation of the extinction from their ratios. 
In the majority of those cases, the extinction Ay is smaller than 1 although in a couple 
of cases it reaches values of up to 4. Note that we find no correlation between Ay and 
the radio luminosity and the luminosity in the lines. There is then no other way than 
to ignore extinction in our work, if we want to work with a sample with significant size. 
Ignoring extinction will then simply add some dispersion to the properties derived from the 
comparison of optical and radio data, which is not really an issue in our work. 

The black hole masses of the galaxies were estimated from the observed stellar velocity 
dispersion given by the SDSS, er*, using the relation by Tremaine et al. (2002): 

logM BH = 8.13 + 4.02 log(o-*/200 km s _1 ). (1) 

In considerations involving black hole masses we disregard cases with a* < 60kms _1 as well 
as cases where the signal-to-noise ratio of the SDSS spectrum at 4000 A is smaller than 10, 
to ensure that the estimate of Mbh is not significantly affected by observational errors. 



3. RESULTS 

3.1. Radio vs. optical luminosities 

Figure 1 presents our sample of radio galaxies in the L 1A - Ln a and L 1A - L[olll] planes. 
Radio galaxies of FRI type are represented by filled red (grey in the printed edition) circles, 
FRII types by filled black circles, and the remaining types, i.e. FRI/II, double-double, X- 
shape and one-sided in open blue (grey in the printed edition) circles. As it is apparent from 
the figure, the radio luminosities strongly correlate with the narrow-line luminosities, with a 
Spearman rank correlation coefficient r$ = 0.67 for Ha lines and r$ = 0.62 for [O ill] lineal 

In the left panel of Figure 1 we also labelled the axes in units of Pj and Ltd, where Pj 
is the jet power and L^oi is the AGN bolometric luminosity. They are obtained by using the 
following conversion formulae: 



2 Rather than the widely used Pearson correlation coefficient to measure the strenghts of correlations, 
we prefer to consider the Spearman rank correlation coefficient which is appropriate for all the diagrams 
presented in this study. Indeed, the use of the Pearson correlation coefficient may lead to spurious interpre- 
tations of correlations in the case where both variables are scaled by a common factor (see e.g. Dunlap et al 
1997, Barraclough 2007). 
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L bol = 2.0 x 10 3 L Ha = 7.8 x lO 36 L„ Q [L ]erg s" 1 , (2) 

(Netzer 2009), and 

Pj = 1.6 x 10 18 (//3) 3/2 Li. 4 [W Hz -1 ]erg s_1 > ( 3 ) 

the latter being taken from Willott et al. (1999) with the following modifications: (i) con- 
version from 151 MHz to 1.4 GHz assuming the spectral index a = 0.8 [L v oc ^ -a ); (ii) 
replacement of Pj oc Lf, by Pj oc L u , the latter taken with a normalization factor giving 
equality of both at L L4 = 10 26 W Hz -1 . With this normalization the modified formula leads 
to overestimation of a jet power by a factor 1.4 for Li. 4 = 10 27 W Hz" 1 and underestima- 
tion by a similar factor for L 14 = 10 25 W Hz -1 . These differences are not substantial when 
compared with uncertainties of an original formula expressed via the parameter /, with its 
1 — 20 range. The figure was made adopting f=3. One can see that even for low values of 
/ the jet powers of many objects exceed their bolometric luminosities. 



3.2. Radio luminosities vs. Eddington ratio 

In Figure 2 we plot the radio vs. emission line luminosities normalized by the black 
hole mass. They correlate similarly strongly as their absolute values, with a Spearman rank 
correlation coefficient r$ = 0.77 and r$ = 0.71 for Ha and [O ill] lines, respectively. Noting 
that Lune/MBH provides the proxy for the Eddington ratio defined to be A = L^j Lem and 
using the conversion formulae and ([3]), one can see in the left panel of Figure 2 that 
NLRG cover about 4 decades of A, from A ~ 10~ 4 up to A = 1. The figure shows that our 
sample is dominated by objects with A spanning the range of 10~ 4 — 10 -2 . These AGN are 
optically too weak to be considered as hidden quasars, but nonetheless, being emitters of 
strong and high excitation lines are expected to be powered like quasars by cold accretion 
disks and, if not obscured by torus, would appear to us as BLRGs (Urry & Padovani 1995). 
However, one should note that the division of AGN for RLQs and BLRGs does not have 
any physical grounds, they form continuous 'Eddington ratio sequence' with no signs of an 
accretion mode change (Sikora et al. 2007) i 

The deficiency of AGN with A > 0.01 in our sample confirms the earlier indications of 
rarity of very high accretion rate sources at low redshifts located in massive galaxies (see, 
e.g., Kauffmann et al. 2008), while the presence of several FRI RGs at A > 0.01 is consistent 



3 Historically the AGN division to RLQs and BLRGs was related to the stellar vs. fuzzy optical appearance 
of host galaxies, presently it is usually related to the specific value of the AGN absolute optical magnitude, 
e.g. Mb = —23.0 as in the Veron-Cetty & Veron catalogs (1993). 
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with a direct finding by Heywood et al. (2007) that radio morphologies of type FRI - which 
are usually associated with low luminosity RGs - do happen in quasars as well. 

3.3. Radio-loudness vs. Eddington ratio 

In Figure 3 we plot the dependence of L 1A /Li ine on L[ ine /M BH , where L 1A /Li ine can be 
considered to be the proxy of the radio-loudness defined as the radio-to-optical flux ratio. 
Our results show a significant negative correlation of radio-loudness with the Eddington 
ratio, with a Spearman rank correlation coefficient r$ = —0.54 and r$ = —0.63 when using 
the Ha and [O ill] lines, respectively. Such an anti-correlation was discovered previously by 
Ho (2002) for radio quiet AGN and by Sikora et al. (2007) for radio-loud AGN, however, 
the statistics of their studies was too poor to claim its presence in the sample when limited 
only to strong-emission-line objects. 

In the left panel of Figure 3 we have also indicated the values of Pj/Lboi and A on the 
axes. We can see that most objects with A < 0.01 have jet powers exceeding bolometric 
luminosity of their AGN, those with lowest A's even by a factor larger than 10 (//3) 3 / 2 (see 
Eq. E]). 

3.4. Sizes 

Figure 4 shows the sizes of the radio lobes as a function of the Eddington ratio obtained 
using Ha (left panel) and [O ill] (right panel). To our knowledge, this is the first time that 
such a diagram is shown. There is only a weak correlation, with a Spearman rank correlation 
coefficient r$ = 0.34 and r$ = 0.27 when using the Ha and [O ill] lines, respectively. This 
suggests that the expansion of radio sources is not accompanied by monotonic, long term 
changes of the accretion rate, and that product of expansion velocity multiplied by the 
source life-time does not depend on the Eddington ratio. Furthermore, the fact that the 
source sizes show no dependence with morphological type suggests that radio morphologies 
do not form any evolutionary sequence. However one cannot exclude "switches" between 
different morphologies caused by modulations of the jet power and the jet direction. 

4. DISCUSSION 

Amongst all AGN, the most spectacular from the observational standpoint, yet most 
challenging theoretically - are probably those associated with extended, luminous radio 
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structures. They appear to us as BLRGs and RLQs if oriented with respect to our line of 
sight such that their nuclei are not obscured by dusty tori, and as NLRG otherwise. Their 
radio and optical luminosities imply an efficient energy transport from nuclei to radio-lobes 
via narrow relativistic jets at a rate very often exceeding the bolometric luminosities of their 
host nuclei (Rawlings & Saunders 1991; Ghisellini et al. 2010; Fernandes et al. 2011; Punsly 
2011; and Section 3 in this paper). 

As it was demonstrated by Tchekhovskoy et al. (2011) and McKinney et al. (2012), such 
powerful jets can be produced in the scenario which involves magnetically arrested/choked 
accretion flows. Such an accretion mode can well take place in the innermost portions of a 
disk, when the amount of the net magnetic flux, $, amassed in the central region is larger 
than the maximal flux which can be imparted on that region around the black hole by the 
ram pressure of accreting plasma, 

^BH = (t>B H R g {cM) l/ \ (4) 

where M is the accretion rate, R g = GMbh/c 2 , and 4>bh is the dimensionless factor called 
by Tchekhovskoy et al. (2011) the "dimensionless magnetic flux." The value of 4>bh depends 
on the details of the model and according to the numerical simulations by McKinney et al. 
(2012) is typically on the order of 50. @ Jets powered by rotating black holes threaded by such 
magnetic flux appear to gain kinetic luminosities (Blandford & Znajek 1977; Tchekhovskoy 
et al. 2010) 

Pj ~ 4.0 x 10" 4 - <$> 2 BH tt% H f(n BH ) ~ 10 (<W50) 2 x 2 J(x a )Mc 2 , (5) 
c 

where 

x a = R g n BH /c = [2(1 + Vl-a 2 )]- 1 a (6) 

f(x a ) ~ 1 + lAx 2 a - 9.2a; 4 , (7) 

and a = J BH j 'Jsu^ax = cJbh/GM bh is the dimensionless angular momentum of a BH, 
commonly named 'spin'. For maximal black hole spins, a ~ 1 (— > x a ~ 1/2), 

P hmax ^ 1.9(<W50) 2 Mc 2 = 19(<W50) 2 W(e/0.1) • (8) 

where e = Lf, i/(Mc 2 ) is the radiation efficiency of an accretion disk. The jet power may also 
contain a contribution from the accretion flow. However, as it was shown by McKinney et 
al. (2012) this contribution is never dominant and therefore will be ignored in our further 



4 Note that dimensionless magnetic flux defined by McKinney et al. (2012) and denoted by T bh is lower 
by a factor 5, and is found to be typically of the order of 10. 
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discussion. As it can be verified using Eqs. (j^J) and (jSJ) and Figures 1 and 2, the model 
predicts domination of jet powers over AGN bolometric luminosities even for moderate spins. 

As the next step, we will investigate whether - and how - such a model can explain the 
anti-correlation of radio-loudness with Eddington ratio shown in Fig. 3. We have 

K = L lA /L Une oc Pj/L bol = Pj/(eMc 2 ) = Vj /e (9) 

where r]j = Pj/Mc 2 is often referred to as the jet production efficiency. Hence the negative 
correlation of 1Z with A may potentially arise from the respective dependencies of rjj and/or 
e on A. According to Eq. fl5]), the dependence of rjj on A can eventually result from the 
dependence of the spin on A. However that would require a negative correlation of spins 
with A, which is the opposite of what one might expect by noting that the black holes are 
span up more efficiently for larger rather than smaller accretion rates. Hence we are left with 
the option that the disk accretion efficiency decreases with the decrease of the Eddington 
ratio, i.e. that e correlates with A. 

The correlation can be explained if assuming that all cold accretion episodes start with a 
similar total net magnetic flux <3>, which is sufficiently large to exceed &bh for any accretion 
rate (see Eq. H]). In such a case, the dynamical dominance of magnetic field over the accreting 
plasma extends up to R m 2> Ri n , where R in is the inner edge of a standard accretion disk not 
affected by magnetic fields, and R m , sometimes referred to as the magnetospheric radius, is 
the size of the region within which the magnetic flux $ is enclosed. Since below that radius 
the accretion proceeds via interchange instabilities, no significant optical-UV radiation is 
expected to be produced within this region. Hence, being 'truncated' at R m , the accretion 
disk is expected to have radiation efficiency e ~ R g /R m . Noting that for a given black hole 
mass R m /R g oc ($) 4//3 m~ 2 / 3 (Narayan et al. 2003) where m = Mc 2 / L EDD = A/e, one can 
find that e oc A 2 / 5 and therefore that 

K oc rjj/e oc A~ 2/5 . (10) 

Some dispersion is expected to be imposed on this relation by distributions of black hole 
masses and magnetic fluxes (unless $ oc Mg^). 

The large ratios R m /R g for AGN with low values of A are consistent with detailed 
observations of some individual BLRGs which strongly suggest existence of the truncation 
radius in their accretion disks (Eracleous et al. 2000; Grandi & Palumbo 2007; Sambruna 
et al. 2009; Tazaki et al. 2010; Cowperthwaite & Reynolds 2012). In those papers, the disk 
truncation radii inferred from observations were interpreted as an effect of obscuration of 
the central region by Thomson thick corona, or as a transition to the advection dominated 
accretion flows. Since it predicts large truncation radii as determined by R m , the model can 



explain the cases of objects with Pj/Lboi > 10 (see Punsly 2011 and refs. therein) without 
the necessity to postulate the jet production efficiency significantly greater than unity@ 



The scenario presented above needs to be verified by demonstrating how such large net 
magnetic flux can be assembled in the central regions of AGN. This problem was raised by 
Lubov et al. (1994), who showed that it is possible only if the magnetic Prandtl number 
is < H/R, where H is the height of the disk at a distance R from the black hole center. 
Since the Prandtl number predicted in fluids with isotropic turbulence is expected to be on 
the order of unity (Parker 1979), the above condition is rather difficult to satisfy in cold, 
geometrically thin disks (Livio et al. 1999; Cuo 2011). This problem can be alleviated, if one 
considers a possibility of an accumulation of a large magnetic flux in the central regions of 
AGN at the onset, or even prior to the period of the cold accretion phase. Such an accumu- 
lation could occur by dragging magnetic fields by geometrically thick, advection dominated 
accretion flows, such as those predicted for super-Eddington accretion rates (Jaroszyhski et 
al. 1980; Beloborodov 1998; Abramowicz 2005) and also for very sub-Eddington accretion 
rates (Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1994). Since our studies of source sizes 
do not indicate any evolutionary trends of expansion of sources with decreasing Eddington 
ratio (Fig. 4), the second option seems to be favored. This is also supported by probabilistic 
arguments. For a given frequency of the cold accretion events, which in turn are very likely 
to be triggered by mergers of giant ellipticals with less massive, cold gas - rich galaxies (see 
Ramos Almeida et al. 2012 and refs. therein), the fraction of those events accompanied by 
production of powerful jets is predicted to correspond to a reasonably significant fraction 
(~ 3 — 10%) of giant ellipticals to be in radio-active states (Sadler et al. 1989; Donoso et 
al. 2009; van Velzen et al. 2012). These ellipticals are most likely powered by the Bondi 
accretion of hot interstellar gas (Burns 1990; Hardcastle et al. 2007; Tasse et al. 2008; Dunn 
et al. 2010; Werner et al. 2012). For typical intensities of interstellar magnetic fields in 
such galaxies, ~ 10/iG, and their coherence scales, ~ 100 parsecs (Moss & Shukurov 1996; 
Mathews & Brighenti 2003 and refs. therein), the accumulation magnetic fluxes in the cen- 
tral regions of such objects needed to exceed $m for any accretion rate can occur within a 
small fraction of the Hubble time. 



Formally r/j can exceed unity, which simply would mean that the rate of extraction of the black hole's 
rotational energy is larger than the rate of the energy inflow (Tchekovskoy et al. (2011). However as most 
recent simulations of McKinncy et al. (2012) indicate, it is rather difficult to achieve such a solution. 
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5. CONCLUSIONS 

Main results of our studies of radio and spectroscopic optical properties of NLRGs at 
z < 0.4 can be summarized as follows: 

• Radio luminosities are found to correlate strongly with narrow line luminosities. When 
converted to the jet powers and AGN bolometric luminosities, they indicate that the jet 
kinetic energy often exceeds the radiative output of accretion flows; 

• NLRGs cover about 4 decades of the Eddington ratio, A, with most of them having A < 0.01. 
This indicates that their parent population is dominated by BLRGs. Together with RLQ 
they form a continuous Eddington ratio sequence; 

• The 'radio- loudness', TZ = Li^/Lu ne , shows a strong negative correlation with A. This 
would imply that the jet production efficiency is the largest at lowest values of A, provided 
disk radiation efficiency is independent on A. 

• The lack of any signatures of correlation or anti-correlation of radio source sizes with 
Eddington ratio indicates the lack of any significant monotonic migration of objects (to 
lower or larger Eddington ratios); 

• A promising scenario which can explain energetics of jets in powerful radio sources and 
observed radio vs. optical luminosity correlations is the one involving the magnetically 
arrested/choked accretion flows. Such flows may support sufficiently large magnetic fluxes to 
power jets with Pj > L^i, while a truncation of accretion disks by a 'poloidal magnetosphere' 
can relax requirements of having rjj > 1. 

• Our results suggest a connection of the cold accretion phase following a lower accretion 
rate, hot accretion phase taking place in extragalactic radio sources. Such a two-phase 
scenario can overcome the difficulty of accumulating large magnetic fluxes by geometrically 
thin accretion disks. Without such a pre-phase, the cold accretion events would not be 
accompanied by production of powerful jets. 
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NASA JPL) no. RSA 1433865. GS and DKW acknowledge financial support from the Eu- 
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Fig. 1. — The radio luminosities L function of the optical line luminosities Lh« (left 

panel) and £[oiii] (right panel). FRI types are represented by filled red circles (grey in the 
printed edition), FRII types by filled black circles, and the remaining types, i.e. FRI/II, 
double-double, X-shape and one-sided in open blue circles (grey in the printed edition). In 
the left panel, the values of the bolometric luminosity and of the radio luminosity Pj in ergs 
s _1 calculated for / = 3 are also indicated. 
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Fig. 2. — The radio luminosities normalized by the black hole masses as a function of optical 
emission line luminosities (left panel: Lho, right panel: £[oiii]) normalized by the black hole 
masses. The symbols have the same meaning as in Fig. 1. In the left panel are also indicated 
the values of the parameters log A and log Pj/L edd for / = 3. 
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Fig. 3. — L 14 /L H q as a function of Ln a /M BH (left panel) and -tq^/I^oiii] as a function of 
-^[oiiii/Mbh (right panel). The symbols have the same meaning as in Fig. 1. In the left 
panel are also indicated the values of the parameters log A and log Pj/L^i for / = 3. 
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Fig. 4. — The projected sizes of the radio lobes as a function of the Eddington ratio, as 
measured using Ha (left panel) and [O ill] (right panel). The symbols have the same meaning 
as in Fig. 1. 



Table 1. Radio and optical properties of the sample of radio galaxies 1 
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1 Table 1 is published in its entirety in the electronic edition of ApJ. A portion is shown here for guidance regarding its form and content. 



